Draft version December 16. 2009 

Preprint typeset using I^TgX style emulateapj v. 08/22/09 



ON THE KENNICUTT-SCHMIDT RELATION OF LOW-METALLICITY HIGH-REDSHIFT GALAXIES 

NiCKOLAY Y. GNEDIn'-' ' AND ANDREY V. KRAVTSOV" ''''' 
Draft version December 16, 2009 



0^ 
O 
O 
(N 

o 

CD 

Q 



O 

U 

6 



ABSTRACT 

We present results of self-consistent, high-resolution cosmological simulations of galaxy formation at z '--^ 
3. The simulations employ recently developed recipe for star formation based on the local abundance of 
molecular hydrogen, which is tracked self-consistently during the course of simulation. The phenomenological 
H2 formation model accounts for the effects of dissociating UV radiation of stars in each galaxy, as well as 
self-shielding and shielding of H2 by dust, and therefore allows us to explore effects of lower metallicities 
and higher UV fluxes prevalent in high redshift galaxies on their star formation. We compare stellar masses, 
metallicities, and star formation rates of the simulated galaxies to available observations of the Lyman Break 
Galaxies (LBGs) and find a reasonable agreement. We find that the Kennicutt-Schmidt (KS) relation exhibited 
by our simulated galaxies at z « 3 is substantially steeper and has a lower amplitude than the z = relation at 
Sh ^100 M0pc"^. The predicted relation, however, is consistent with existing observational constraints for 
the z « 3 Damped Lyman a (DLA) and LBGs. Our tests show that the main reason for the difference from 
the local KS relation is lower metallicity of the ISM in high redshift galaxies. We discuss several implications 
of the metallicity-dependence of the KS relation for galaxy evolution and interpretation of observations. In 
particular, we show that the observed size of high-redshift exponential disks depends sensitively on their KS 
relation. Our results also suggest that significantly reduced star formation efficiency at low gas surface densities 
can lead to strong suppression of star formation in low-mass high-redshift galaxies and long gas consumption 
time scales over most of the disk in large galaxies. The longer gas consumption time scales could make disks 
more resilient to major and minor mergers and could help explain the prevalence of the thin stellar disks in the 
local universe. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: formation - stars:formation - methods: 
numerical 
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1. INTRODUCTION 

The hierarchi cal Cold Dark Matter (CDM) structure forma- 
tion paradigm ( White & Reeslll978t iBlumenthal et alj|1984 
has proved to be remarkably successful both in predicting and 
explaining a variety of observational data from the detailed 
properties of the temp erature an isotropics of the Cosmic Mi- 
crowave Background jDunklev et al.ll2009l) to the linear (e.g. . 



Tegmark et al. 120061) and nonlinear (e.g., Springel et al|2006 
Conrov et all l2006riKravts6vl llOOa iWetzel & Whitd I200S 



clustering of galaxies. 

Galaxy formation modeling within the CDM frame- 
work has produced a general picture of how galaxy- 
sized objects collapse out of primo rdial density fluctuations 
(e.g.,[S teinmetz & Navarro 2002; S pringel & Hernquisll20"03l: 
iMaver et al.. .2008;) . Nevertheless, so me key aspects of the 
picture are still b eing developed ( Birn boim & Dekel 120031 : 
iKeres et al]l2005 ^ and accurate predictions of even the basic 
properties of galaxies, such as luminosities, col ors, metallic- 
ities, and gas content remain a challenge (e.g., iMaver et alj 
120081) . These properties depend sensitively on the details of 
star formation. The root of the challenge is therefore our lack 
of full understanding of how gas is converted into stars under 
different conditions during different stages of galaxy evolu- 
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tion. 

Empirical studies of the Kennicutt-Schmidt (KS) relation 
([Schmidt 1959; Kennicutt 1989, 1998) between the total sur- 
face density of the neutral gas Sh, or its specific phase (atomic 
or molecular), and surface density of star formation rate Ssfr 
established that the KS relation can be approximated by a 
power-law, Ssfr oc S^, with n « 1.4 (or, more generally, be- 
tween 1 and 2) at high surface densities, when Sh is dom- 
inated by the molecular gas, but steepens considerably at 



at Sh < 10M(r,/ pc2 (iK ennicuttI Il998t iMartin & Kennicutt 



1 2001 1 : [Wong & Blitz 2002 at iHever et al.l 2004: Boissier et al 
hoOT-. i Bigiel et aL,200& iRovchowdhurv et al.i,2009i: .Bolatto 
2009). These studies, however, have been limited to samples 
of nearby galaxies with a relatively narrow range of metallic- 
ities (« 0.5- lZ(n) and mo rphologies (thin gaseous disks; see 
however fRovchowdhurv et al. 2009). 

Observational picture remains sketchy at higher redshifts 
and lower metallicities. On one hand, recent studies show that 
the power law KS relation between star formation and molec- 
ular gas surface density hol ds in actively s tar forming galax- 
ies at higher redshifts (e.g.. iBouche et al.H2007l : iBaker et^ 
l2004al) . On the other hand, recent searches for the UV emis- 
sion as sociated with young stars in damped Lya systems 
(DLA; I W olfe & Cherj j2006b and in Call absorbers at low 
redshifts (Wild et al. 2007) indicate that Esfr at the gas sur- 
face densities characteristic for these systems (A^hi ~ 10^' - 
10^^ cm"^ or Eh ~ lO-lOOMo/pc^) is an order of magnitude 
lower than implied by the local KS relation. A similar result 
is indicated by recent studies of the Lyman Break Galaxies 
dkafelski et al. 2009; Rafelski 2009). 

Given that the typical metallicity of the DLA systems is 
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considerably lower than solar^ dFox et alj 120091 and refer- 
encse therein), this may indicate that the Esfr - Sh relation is 
significantly lower at these surface densities in systems of low 
metallicity. The first evidence of this in the nearby galaxies is 
that the KS relation in the SMC appears to be significantly 
steeper than « w 1 .4 at Sh < SOMq/ pc" (Bolatto 200 9). 

Overall, the trends described above can be expected. Ob- 
servations show that star formation in galaxies correlate s 
strongly with the molecular gas (e.g.. lWong & Blit3l2002bl) . 
At the same time, SFR has only weak or no correla- 
tion with the d ensity of atomic gas (Wong^^Blitz 2002b; 
iKennicutt elaH i2007: Bigiel et al. 2008D. We can thus ex- 
pect that Esfr relation with Eh = Eh, + Ehi will depend on 
the molecular fraction of the gas /h, = Ehj/Eh- The molec- 
ular fraction in observed galaxies is controlled by the gas 
density, FUV flux, and pressure (e.g.. Blitz & Rosolowsk^ 
2006), con sistent with theoret ical expectations (Elmegreeji 
1993: Elmegreen & Parravanol il994 ; Ro bertson & Kravtsov 
20081) . but is also expected to be sensitive to the dust con- 
tent (and thus metallicity) of the gas because dust plays an 
important role both in shielding molecular ga s from UV radi- 
ation and in catalyzing production of H2 (e.g.. Sta bler & Pallal 
W005). Mor cover, the efficiency with which a galaxy is able to 
build high-density molecular regions should depend on metal- 
licity, because the compression of gas in the radiative shocks, 
arising in the highly turbulent interstellar medium (ISM) of 
gaseous disks, should depend on the metaUicity-dependent 
cooling rate. 

These considerations suggest that in order to make progress 
in our modeling of star formation in galaxies, we need simula- 
tions with both the high spatial and mass resolution to resolve 
the density structure and thermal state of different regions of 
the ISM in galaxies and a model for formation of molecular 
hydrogen. Several theoretical models incorporating some 
of these elements and potentially capable of predicting 
dependence of star formation relati on on ISM proper- 
ties have bee n devel oped recently (Krumholz & McKed 
2005; Pelupess vetaTI [2 006; Robertson & Kravtsov 2008; 
Tasker & Brvan '2008; Schave & Dafla Vecchia 2008; 
&iedinetal. 2009; Tasker & Tan 2009; Krumholz et al. 
2009btlPelupessy & Papadopou los 2009). These models have 
shown that thermodynamics of the ISM can control the shape 
of flie density PDF of th e ISM gas (Wada & Norman 2001, 
2007; 'Wadaetal.' "^OO?, 'Kravtsov) i2003i: ,Tasker & Brvam 
2008; Robertson & Kravtsov 2008) which, in turn, can affect 
the KS relation. The general prediction of the models that 
incorporate the physics of atomic-to-molecular transition at 
the boundaries of molecular clouds is that the density of this 
transition, controlled primarily by dust abundance (i.e., metal- 
licity of the gas), is a ke y factor affecting the amphtude and 
slope of the KS relation (Gnedin et al. 2009; Krumhol z et al.l 
|2009b). The models of Robertson & Kravtsov (2008!) have 
also highlighted importance of the interstellar radiation 
field for contr olling the abundan ce of low-density, diffuse 
molecular gas (lEknegreenl 1 1 9931) . as wefl as the slope and 
amplitude of the large-scale KS relation and metallicity of 
the gas. So far, however, none of the models have considered 
ISM of galaxies arising and evolving self-consistently in 
cosmological context. 

In this paper, we present results of self-consistent, 

^ We refer to "solar metallicity" as a typical gas metallicity in the solar 
neighborhood, Zq = 0.019 by mass, rather than the metallicity of the Sun, 
which has been recently revised downward. 



high-resolution cosmological simulations of galaxy for- 
mation, which em ploy the metaUicity-depen d ent m odel 
of H? presented in Gne din, Tassis. & Kravtsov! (12009 1) and 
iGnedin & Kravtsov (2010). In the latter paper we also study 
sensitivity of global star formation to the full range of metal- 
licities and UV fluxes. The focus of this paper is to explore 
the predictions of our H2-based star formation recipe for the 
KS relation in the low-metallicity, high UV-flux environments 
of high-redshift galaxies and compare these predictions with 
the existing observational constraints from DLA and Lyman 
Break Galaxies (LBGs). 

2. SIMULATIONS AND STAR FORMATION MODEL 

The physical ingredients and computational setup for the 
simulation we use in this paper are described in lGnedin et al.l 
( l2009i) . As a brief mmmary, the Adaptive R efineme nt Tree 
(ART ) code (iKravtsov 1999; K ravtsov et al.l2 002; Rud d et all 
l2008h uses adaptive mesh refinement in both the gas dynamics 
and gravity calculations to achieve high dynamic range in spa- 
tial scale. The high-mass resolution is achieved by using the 
standard "zoom-in" initial conditions with particles of smaller 
mass resolving Lagrangian region of a system of interest. 

The specific initial conditions setup we use in this paper 
is a Lagrangian region corresponding to five virial radii of a 
galactic system, which evolves to approximately Milky Way 
mass (M « 10'' Mq) atz = 0, with the mass resolution of 1.3 x 
10'' M0 in dark matter, 2.2 x 10^ M© in baryons, and with 
the spatial resolution of 65 pc x [4/(1 +z)] (in physical units) 
within the region. This Lagrangian region is embedded into 
a cubic volume of 6/;"' comoving Mpc on a side, which is 
followed with a coarse 64^* grid outside the Lagrangian region 
with periodic boundary conditions. 

Our fiducial simulation includes star formation and super- 
nova enrichment and thermal energy feedback, and follows 
self-consistently the 3D radiative transfer of UV radiation 
from in dividual stellar particl es using the OTVET approxi- 
mation ( iGnedin & Abelll2001h . The simulation incorporates 
non-equilibrium chemical network of hydrogen and helium 
and non-equilibrium, metaUicity-dependent cooling and heat- 
ing rates, which make use of the local abundance of atomic, 
molecular, and ionic species and UV intensity, followed self- 
consistently during the course of the simulation. 

A phenomenological model of molecular hydrogen forma- 
tion is used to identify the locations of molecular clouds in 
the simulated galaxies, as described in Gnedin et al.' ('2009!) 
and, in more detail in Gnedin & Kravtsov (2010). The model 
incorporates both self-shielding of molecular hydrogen from 
the UV radiation and the shielding by dust and takes into 
account catalyzing effect of dust on H2 formation. The lo- 
cal abundance of dust in our model is assumed to linearly 
scale with gas metallicity, which is co nsistent with obser- 
vations of Qie MW LMC , and SMC dGnedin et al.l 120091: 
IGnedin & Kravtsovl 1201 Oh . Gas metallicity in the simula- 
tions is modeled by including metal enrichment by super- 
novae (both type II and la) assuming standard yields, and 
by stars via wind mass loss, with the total amount of heavy 
elements released into the ISM determined by the assumed 
Miller-Scalo IMF of the stars (see K ravtsov..20 03. for details). 
Once released into the gas the heavy elements are advected by 
the code in the same way as the gas. 

We calibrate our H2 model to fit the observed atomic and 
molecular gas fraction as a function of column density in the 
different metallicity environments of the Milky Way, LMC, 
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and SMC and other nea rby galaxies dGnedin et alj l2009t : 
iGnedin & Kravtsovll2010l) . The model reproduces the metal- 
licity dependence of the column density at which a sharp tran- 
sition from the atomic to fully molecular gas occurs. 

The star formatio n recipe adopted in this work is similar 
to the recipe 2 of Gnedin et alJ (1 20091) . albeit with somewhat 
different parameters. Specifically, stellar particles are formed 
in cells that have mass fraction of molecular hydrogen higher 
than /h, > 0. 1 according to the rate: 



£ff 

P* = —Piil 



(1) 



where is the density of molecular hydrogen, Tsf = 
min(Tn,ax, Tf{), Tff = \Jl>n jlilGp is the free-fall time for a uni- 
form sphere, T^ax = 6.65 x 10^ years is the free-fall time in 
the gas with the total hydrogen number density of 50 cm"^, 
and £fF is star formation e f ficiency per free fall tim e (e.g., 
iKrumholz & McKed l200l |McKee & Ostriker 2007). We 
adopt the efficiency of Eff = 0.005 consistent with observa- 
tional constraints for the average efficiency of star forma- 
tion in the mole cular clouds of the Mi l ky Way and other 
nearby galaxies (IKrumholz & TanI l2007t iMcKee & Ostrikeil 
l2007ft . The Tsf we adopt assumes that in low density cells, 
in which molecular fraction /h, is below unity, star formation 
proceeds mainly in unresolved molecular clouds on subgrid 
scales. This assumption then also motivates setting the maxi- 
mum free fall time to T^ax corresponding to the number den- 
sity of 50 cm"^ typical average density of molecular clouds. 
The /h, < 1 in these cells then can be viewed as reflecting the 
fraction of the total gas in such star forming molecular clouds, 
which themselves have /hj = 1, rather than incomplete con- 
version of the atomic gas into the molecular form inside the 
clouds. 

Note that the empirical calibration of both the H2 forma- 
tion model calibration and the efficiency of star formation is 
done on small scales of individual molecular clouds. The re- 
sults presented in this paper, however, are on global kiloparsec 
scales and the model was not adjusted in any way to produce 
these results. 

3. STAR FORMATION IN HIGH REDSHIFT GALAXIES 

3.1. Global properties of simulated galaxies 

Before we address the star formation relations, we need to 
consider the bulk properties of simulated galaxies to make 
sure they can be plausible counterparts of observed LBGs and 
DLAs. In addition to the main Milky Way progenitor, the 
fully resolved Lagrangian region of the simulation volume 
contains a number of smaller mass galaxies. Figure [T| shows 
the stellar masses of the simulated galaxies versus the total 
mass of their parent halos, defined as a mass within a spher- 
ical radius enclosing overdensity of 180 with respect to the 
mean density, at redshifts 2.5 < z < 3.5.^ The figure shows 
that simulated galaxies span almost two orders of magnitude 
of total mass (Mtot 10'°-6 x 10" Mq) and more than two 
orders of magnitude of stellar mass (M, -lO**-5xlO"'M0). 
Interestingly, the M» - Mtot correlation has the shape quali- 
tatively similar to that required to explain the difference be- 
tween the predicted halo mass function and observed stel- 
lar mass function at lower redshifts (e.g., Conroy & Wechsler 

' We only include in Fig.[T]and the following figures galaxies whose cen- 
ters are fully resolved to the highest resolution of our simulation (65 pc at 
z = 3, nine levels of refinement), with the sole exception of Fig.|2] 




Mtot (Mo) 

Fig. 1 . — The stellar mass of fully resolved model galaxies versus the total 
mass of their halos. The solid black line is the limit of the complete conver- 
sion of baryons into stars, M, = (f2g/r2M)Mtot- The sets of colored points 
correspond to the three redshifts in the simulation that cove r the range of 
redshifts probed by observed systems of lWolfe & Chenl 120061) . 



120091: iMoster et al. 2009b). This indicates that processes in- 
cluded in our simulations could, at least partly, be responsible 
for the observed shape of stellar mass function. We note, how- 
ever, that normalization of this relation in our simulations is 
considerably higher tha n that indicated by observ ational ana- 
lyzes (e.g., see Fig. 2 in lConroy & W echslei'2009^. especially 
for lower mass systems. This indicates that star formation 
in small-mass halos is still considerably more efficient in our 
simulations than in real galaxies. 

The star formation rates, stellar masses, metallicities, and 
radiation fields of these galaxies are shown in Figures |2] at 
three redshifts in the range 2.5 < z < 3.5. For comparison, we 
also show for the observational estimates of the same prop- 
erties f o r the Lyrnan Brea k Ga laxies from ^annucci et alj 
( I2009bl) . lErb et al.1 ( l2006ah . and iPettini et alj (l2001h . Since 
the volume of our simulation is relatively small, the overlap 
with the observed galaxies in stellar mass is limited. In the 
region of overlap, however, our model galaxies have star for- 
mation rates and stellar masses consistent with those of ob- 
served LBGs. The gas metallicities of our model galaxies are 
consistent with observations within the errorbars^. 

The average metallicities of the simulated galaxies are 
higher than the typical metallicity of DLA systems at z ^ 3 
(Prochaskaet al. 2003; Fynbo et al. 2008). This may be be- 
cause the masses of most DLA systems are lower than the 
range of halo masses resolved in our simulations or because 
DLA systems probe the outer regions of the disks where 
metallicities are systematically lower due to metallicity gra- 
dients. It is also possible that we overestimate star formation 
efficiency and hence metal enrichment in objects of a given 
mass. The trends of star formation relations with metallic- 
ity discussed below can thus be expected to be even stronger 
for the real DLA systems compared to the higher-metallicity 
galaxies in our simulations. 

The radiation field in a galaxy may vary widely between 
different locations and environments. To estimate a "typi- 
cal ISM value", we show the median radiation field and its 
10% to 90% range in the gas with density between 1 cm"-' and 
100 cm"-' in the bottom panel of Fig. |2] The radiation fields 
in the ISM of our model galaxies are significantly above the 

^ We note, however, that metallicities in the simulations carry all the sys- 
tematic uncertainties of the adopted nucleosynthetic yields and of the IMF. 
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Fig. 2. — The gas metallicity (top), the stellar mass (middle), and the radia- 
tion field (bottom) of fully resolved model galaxies versus their star formation 
rate. The error bars on the metallicity and radiation field points show the 10% 
- 90% range of these properties measured in the individual galaxies. The sets 
of colored points correspond to the three redshifts in the simulation that cover 
the range of redshifts probed by observed systems of Wolfe & Chen 1 2006). 
Black filled circles, crosses, and open squares are the observational measure- 
ments for LB G galaxies from Mannucci et al. 1 2009b), Erb et al. 1 2006a), and 
IPettini et alj (2001 ,), respectiv ely. Black squares are measurements for GRB 
host galaxies from lChen et al] (20091) . 



typical UV intensity in the Milky Way (which we round off 
to /mw = 10* photons/ cm^/ s / ster/ eV at A = 1000, consis- 
tent w ith the measurements of lDrainelll978l and lMathis et al.l 
11983b . but in agreement with t he typical radiatio n fields mea- 
sured in high-redshift galaxies jChen et al.ll2009t) . 
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Fig. 3. — The KS relation for the simulated galaxies at z = 3 (red): the 
solid line shows the mean relation, while the hatched band shows the RMS 
scatter around the mean. The total surface density Sh = Ehi + Shj takes 
into account only neutral gas. The dotted and short-dashed lines show the 
mean KS relations for the atomic and molecular hydrogen, respectively. A 
solid blue line shows the mean relation only for gas with metaUicities below 
0. IZq , which more closely corresponds to the range of observed metaUicities 
in the DLAs. The long-dashed line is the best fit relation of Kennicutt 1 199^ 
for z ~ galaxies. The gray shaded area shows the measurements of the 
KS relation in LBGs (Rafelski et al. 2010, in preparation). The solid circle 
at high gas surface denstity shows measurements of SFR and molecular gas 
surface densities for the lensed z fa 3 LBG cB58 (Baker et al. 2004a). Solid 
squares and triangles show observational upper Hmits for the z 3 DLA sys- 
tems from Wolfe & Chen (2006), assuming two geometrical models for the 
atomic gas distribution in the DLA systems (squares: uniform disks; trian- 
gles: isolated clouds). The three different points of each type correspond to 
the three repre sentative choices of the smoothing angular scale in Table 1 of 
IWolfe & ChenI (2006) . Errorbars on the triangles reflect the current uncer- 
tainty in the high-end slope of the DLA column density distribution. 

3.2. The Kennicutt-Schmidt relation at high redshifts 

Figure [3] shows the KS relation between star formation and 
gas surface densities at z = 3 for galaxies formed in our cos- 
mological simulations. The star formation rate is averaged 
over 20 Myr, and the gas and SFR surface densities are mea- 
sured on the scale of 500 pc. We have verified that our re- 
sults are robust to changes in the averaging scale (as long 
as the scale is > 200 pc) and the period of time over which 
the star formation rate is averaged (for time intervals smaller 
than 30 Myr). The figure shows that the predicted KS relation 
for the z = 3 galaxies is s ignificantly steep er than the rela- 
tion for z = galaxies (the iKennicutal 19981 fit to local galax- 
ies is shown by the dashed hne) at Sh ^ 100 M©/ pc^ and 
its amplitude is an order of magnitude lower than the ampli- 
tude of the iKennicuttI ( 119981) fit at Sh < SOMg/pc^. The 
corresponding Ssfr-Eh, and Ssfr-Shi correlations in Fig- 
ure [3] show that the steep Esfr - Sh relation is due to the 
rapidly decreasing ratio of molecular to atomic gas, Sh, /Shi, 
at Sh < lOOMo/pc^. Likewise, the surface density at which 
HI surface density saturates is « 50 Mq/ pc^, considerably 
larger value t han saturation surface density observed for local 
galaxies (e.g.. lWong & Blitzll2002al) . 

The predicted KS relation for z ^ 3 galaxies is consistent 
with the recent measurements of star formation densities in 
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Lyman Break Galaxies by Rafelski et al. (2010, in prepara- 
tion), which are shown as a grey band in Fig. [3] and are com- 
plemented by the measurement of Ssfr and Sh for the lensed 
z « 3 Lyman Break Galaxy cB58 (Baker et al. 2004a) at high 
gas surface density. 

Our simulations are also consistent with the upper lim- 
its on star format ion surface density in the z « 3 DLA sys- 
tems obtained by IWolfe & ChenI (120061) . which are shown in 
Fig. [3] with solid triangles and squares. The three different 
points of each point type correspond to the three represen- 
tative choices of the smoothing angular scale in Table 1 of 
[Wolfe & Chen (2006). Solid squares show the original limits 
of 1 Wolfe & Chen (2006.) . However, these limits were derived 
assuming a specific model of "monolithic" gaseous disks for 
DLA systems. This model predicts a specific slope of -3 
on the high end of the HI column density distribution. This 
value of the slope is only marginally (at 3cr ) consistent with 
the value of -6.4+| g measured by Prochaska & Wolfd (120091) 
from the Data Release 5 (DR5) of the Sloan Digital Sky Sur- 
vey (SD SS), but is close to the v alue of -3.5 measured re- 
cently bv lNoterdaeme et al.l (l2009h . The two recent measure- 
ments indicate that there remains a substantial uncertainty in 
the high-A^ slope of the DLA column density distribution. 

Given the uncertainty in the high-end slope, we also con- 
sider an alternative model for DLAs, in which DLA absorp- 
tion arises in a single spherical cloud. The clouds may con- 
tain a mixture of atomic and molecular gas and may have an 
arbitrary spatial distribution (including the disk-like), as long 
as a probability of intersecting two clouds by a single line of 
sight is small. This picture may correspond, for example, to 
the DLAs arising in the extended HI envelopes of molecular 
clouds. In this model any value for the high-end slope of the 
column density distribution is allowed. These two geometric 
DLA models bracket the plausible range of possible geome- 
tries a nd, thus, account for t he implicit horizontal uncertainty 
of the IWolfe & ChenI (|2006l) upper limits. 

These limits, recomputed for the single cloud model, are 
shown in Fig.|3]as downward triangles, with the width of the 
errorbars on the triangles reflecting 3cr errors on the observed 
high-end slope from Prochaska & Wolfe (2009). The error 
bars also encompass the Noterdaeme et al.l (120091) value and 
the value of -3 for the monolithic disk model. The error bars 
of the triangles may serve, therefore, as a reasonable estimate 
of the horizonthal uncertainty of the limits due to the current 
uncertainty in the hig h-end DLA column de nsity distribution. 
We have converted the lWolfe & Che n (2006) upper limits into 
Units in the Ssfr - plane by solving their equations 3 and 
6 (using the values from Table 1) for normalization of the KS 
relation, K, and effective DLA column density, A^, and con- 
verting into Eh- 

This conversion itself relies on a subtle assumption that re- 
quires clarification. DLA lines of sight probe column density 
averaged on transverse scale of ^ 1 pc, while Sh in Fig.[3]is 
averaged on the scale of 500 pc. The relation between the gas 
column densities would of course be direct if the ISM is uni- 
form over this range of scales, but the ISM density both in real 
and simulated galaxies is very far from uniform and exhibits 
densities ranging by several orders of magnitude from the av- 
erage density regions (n ^ 0.1 - 1 cm""*) to dense molecular 
clouds (n > 100 cm"^). It may thus seem like there may not be 
any connection at all between the column densities measured 
by the DLAs and average Eh, which is a mass-weighted av- 
erage measured on '--^kpc scale. For example, imagine a sim- 



ple toy ISM model, in which most of the mass is in dense, 
massive clouds with a covering fraction so low that it is very 
improbable to intersect them with DLA los. In this case, area- 
weighted surface density can be arbitrarily low, while Eh av- 
eraged on 500 pc scale can be arbitrarily high. 

However, there are at least two arguments that correlation 
between and Eh should be tight even for non-uniform ISM. 
First, the DLA column density is an average over tens and 
hundreds of pc in the los direction to quasar Second, av- 
erage effective A^ estimated for an ensemble of DLAs corre- 
sponds to the area-weighted average surface density of ISM 
in galaxies probed by the DLA lines of sight. For a real- 
istic probability density distribution (PDF) of ISM the area- 
weighted mean surface density should be quite close to the 
mass-weighted average. For example, for the log-normal pro- 
jected density PDF, which is a sensible approximation to the 
PDF of real and model galaxies (e.g., Elmegreen 2002), the 
area- and mass-weighted surface densities are directly related 
(see, e.g.. lOstriker et al] 1200 lb . We note, however, that the 
exact form of the PDF for real galaxies is yet known and so 
the factor in convertin g A^ into Eh is uncerta in. Our conver- 
sion and conversion in IWolfe & ChenI (12006b . assume a tight 
one-to-one relation between A^ and Eh. 

Figure [3] shows that the KS relation in our z = 3 galaxies 
is gene rally consistent with the upper limits of IWolfe & ChenI 
( 120061) . especially if we recall that metallicities of these galax- 
ies are somewhat higher than the typical metallicities of 
DLAs. If we only consider the gas with Z < O.IZ0, typi- 
cal for the DLA systems, our predictions are fully consistent 
with the Wolfe & Chen (2006) upper limits. 

What is the primary cause of the difference between z = 
KS relation and the z = 3 relation indicated by observations 
and our simulations? The high-redshift galaxies differ from 
their low redshift counterparts by their systematically lower 
metallicities and higher interstellar UV fluxes.** The low 
metallicities are due to generally lower stellar masses and 
higher gas fractions of high-redshift galaxies, which are in 
the earlier stages of chemical evolution compared to their lo- 
cal counterparts. The higher UV fluxes are due to the con- 
siderably higher star formation rates, smaller sizes, and lower 
abundance of the UV absorbing dust in high redshift systems. 

In principle, both lower metallicity and higher UV flux 
can affect star formation in galaxies. Lower metallic- 
ity and dust abundance makes it more difficult t o form 
dense molecular regions in w h ich star s form^ (e.g., Schavd 
'2004'; Ki-umholz et al. 2001 l2009al: iGnedin et al.l 120091: 
Gnedin & Kravtsov 2010), whi le high UV flux can dissociate 
such regions (e.g., |E lmegreen! ll993l : iRobertson & Kravtsovl 
l2008t IGnedin & Kravtsov.201(f ! 

To gauge the relative importance of the metallicity and in- 
terstellar UV flux we have rerun our simulation from z = 4 to 
z = 2.5 (« 1 Gyr) resetting the dust-to-gas ratio used in the 
calculations of H2 abundance and UV field to fixed multiples 
of the Milky Way values. We present the results of such test 
simulations for a full grid of metallicities and UV intensities 
in the companion paper dGnedin & Kravtsovll2010l) . Here we 

There are of course further differences between local and high redshift 
galaxies, such as geometry of the gas distribution, amount of ISM turbulence, 
etc. These, however, appeal' to have relatively minor effect on the KS relation. 

' For example, the transition from atomic to fully molecular gas occurs 
at A'h ~ 3 X 10^" cm"^ in the Milky way, but at Wh 3 x 10-' cm"^ in 
the factor of 5-10 lo wer metallicity ISM of the SMC ( Browning et al. 200|; 
IGillmon et alJ20(M this difference is partly due to a different dust extinction 
law). 



6 




1 10 102 



^H. Shi. (M^pc^) 

Fig. 4. — The KS relation for the simulated galaxies at z = 3, but with 
dust-to-gas ratio and UV flux fixed to their solar neighborhood values. The 
fines and cross hatched band have the same meaning as in Fig.[3]and the total 
surface density Eh takes into account only neutral gas. The gray shaded area 
shows the KS relation for the local dwarf and normal spiral galaxies measured 
by the THINGS project iBigiel et al. 2008). 

present only results for the rerun in which dust-to-gas abun- 
dance and UV flux were set to their Milky Way values. Fig- 
urelHshows the KS relation for this test re-simulation at z = 3. 
As before, the star formation rate in the simulations is aver- 
aged over 20 Myr and the gas and SFR surface densities are 
averaged on the scale of 500 pc. The figure shows that in the 
case of the MW metallicity and UV flux the predicted KS rela- 
tion i^_consktentwathJhe^ for the z = galax- 
ies ("Kennicutt" 1 998'-. 'Bigi el et al.ll2008l) . In particular, the pre- 
dicted KS relation matches the approximately linear relation 
relation Ssfr oc S]|i° at Eh > lOMo/pc^ (with indication of 
steepening to Esfr oc E^'^^' '' at Eh > IO^Mq/pc^) and the 
significant steepening of the relation at Eh < 10 MqPC"^. 
Simulations also reproduce the observed scatter around the 
mean relation reasonably well. Finally, note that the atomic 
hydrogen in this case saturates at surface density of Ehi ~ 
10 MqPC"^, which is consistent with observations of nearby 
galaxies (e.g.. lWong & Blitz 2002a; Bigiel et al. 2008) and a 
factor of five smaller than the corresponding value for the self- 
consistent metallicities and UV fluxes of the simulated z = 3 
galaxies in Fig. [3] 

The fact that the local KS relation is reproduced in simula- 
tions of high redshift disks with higher metallicity and lower 
UV flux indicates that the differences in the ISM structure be- 
tween high and low redshift galaxies play a relatively minor 
role in shaping this relation. Although the density PDF in 
gaseous disks ca n be important in shaping the slope of the KS 
relatio n (iKravtso v 2003; Tassis 2007; Robert son & Kravtsovl 
12008*). our test shows that the main factors behind the differ- 
ence between z = and z = 3 KS relations are gas metallicity 
and UV flux in the ISM. Furtherm ore, additional experiments , 
presented in the companion paper ( Gnedin & Kravtsovl2010l) . 
show that the difference in metallicity is the primary reason 
behind the differences in the KS relations shown in figures |3] 
andH 



The UV flux does affect star formation in the individual 
patches of the ISM significantly, but it also affects the surface 
density of atomic and molecular gas. The changes of Esfr 
and Eh = Ehi + Eh, affected by the change in the UV flux are 
such that a given patch of ISM is moving approximately along 
the KS relation for a varying FUV flux. The average KS rela- 
tion is therefore less sensitive to the changes in UV flux than 
to the changes in the metallicity of the same magnitude, even 
though star formation in individual ISM regions and galaxies 
is strongly affected. The change of metallicities, however, af- 
fects Esfr without affecting the surface density of neutral gas 
significantly, thereby changing the mean KS relation. The ef- 
fect of the UV flux on the KS relation does increase for lower 
metallicity galaxies, as we discuss in the companion paper 

The conclusion from our experiments is therefore that the 
main reason the KS relation of high redshift galaxies is signif- 
icantly different from their low-redshift counterparts is their 
lower metallicity and dust abundance, while the higher FUV 
flux and differences in density PDFs are secondary factors. 

4. DISCUSSION AND CONCLUSIONS 

In this study we use recently developed recipe for star 
formation in simulations of galaxy formation based on the 
local abundance of molecular hydrogen (eq. [U to exam- 
ine the Kennicutt-Schmidt relation expected in high redshift 
systems. The H2 abundance is tracked self-consistently in 
the simulation using a phenomenological model, which in- 
cludes effects of dissociating UV radiati on, as well as self- 
shielding and shielding of H2 by dust (see lGnedin et al.ir2009l; 
[Gnedin & Kravtsov 2010, for details of star formation recipe 
and H2 model). 

Our main result, presented in Fig. |3] is that the KS rela- 
tion in the low-metallicity and high UV flux environments of 
z ~ 3 is substantially steeper and has a lower amplitude than 
that of the z = galaxies at Eh < 100 M0 pc"-. As discussed 
in the previous section, the main reason for the difference is 
lower metallicity and dust-to-gas ratio of high redshift galax- 
ies. While our tests show that the UV flux does affect the star 
formation drastically, it simultaneously affects surface density 
of neutral gas in such a way that the shape and normalization 
of the KS relation are largely unaffected. 

Our predicted KS relation for z ~ 3 galaxies agrees well 
with a recent estimate of Esfr in the Lyman Break Galax- 
ies by Rafelski et al. (2009, in preparation) at low gas surface 
densities and measurements of SFR and gas surface den sity in 
LBG cB58 at high surface densities (Ba ker et af]l2004al) . Our 
results are also consistent with the upper limits on star for- 
mation surface density derived for the subset of z 3 DLA 
systems at gas surface densities of Eh ^ 20-50 M0 pc"^ 
dWolfe & Ch en 2006). These limits are lower by a factor of 
ten compared to the Esfr for the local galaxies at the same 
gas surface densities (Kennicutt 1998) similar to our simula- 
tion results. We note, however, that the typical metallicity of 
the gas shown in Fig.|2]is a few times higher than the typical 
metallicity of the DLA systems (e.g., iProchaska et al.ll2003l; 
iFvnbo et"ani2008 ). Given that effect of decreasing metallicity 
on the KS relation is to steepen it at low surface densities and 
lower its amplitude, we can expect that for systems of lower 
metallicity than those we have considered in this study the KS 
relation should have even steeper slope and lower amplitude, 
as is shown by the solid blue line in Fig. [3] 

A qualitatively similar trend of the KS relation with metal- 
licity was recently predicted by Krumhol z et al.l (12009b') using 
a model of atomic to molecular transition in molecular com- 
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Fig. 5. — The face-on view of the disk of young stai's of the most mas- 
sive galaxy in the simulation at z = 2.5. Only stars younger than 50Myr are 
shown; the colors are indicating the B-V color (bright white con'esponding 
to the bluest, youngest stars). The vertical size of each panel is identical and 
is approximately equal to 10 physical kpc. The top panel shows result of 
the self-consistent simulation (i.e., the simulation for which the KS relations 
shown in Fig. [5), while the bottom panel shows result of the test run, in which 
dust abundance and FUV flux were set to their Milky Way values at ^ = 4 and 
the simulation was run with these tixed values to z = 2.5 (KS relation for this 
run is shown in Fig.|4). Note that star formation in the self-consistent case is 
concentrated towards higher gas density regions, which results in more com- 
pact distiibution of stars in the disk and fewer young stars in the inter-arm 
regions, as shown in Fig.|6] 

plexes based on the lWolfir e et alJ ('2003') semi-analytic model 
of atomic ISM (see also McKee & Krumholz 2009). Their 
model predicts significant steepening of the KS relation be- 
low gas surface density of Sh < lO/cZM0pc~^, where Z 
is the metallicity of the gas in the units of Zq and c is the 
ISM clumping factor. The latter reflects the difference be- 
tween the surface density of the ISM averaged on some scale 
> 100 pc and the surface density of individual giant molec- 
ular complexes on the scale of 100 pc. For averaging scales 
of 500 pc us ed in our calculations the clumping factor should 
be c ~ 2-5 (iKrumholz et al..,2009h) . The model thus predicts 
the steepening of the KS relation at Eh ~ 2-5 M© pc"^ for so- 
lar metallicities and Eh ~ 6- 15M0pc"^ forZ « O.2-O.5Z0 
characteristic for our simulated galaxies. These values are 
considerably lower than the transition surface density of Eh ~ 
50- lOOMopc"^ we derive for our simulated galaxies. It is 
not entirely clear what the source of this difference is. We 
note, however, that the mass fraction of molecular gas for a 
given large-scale surface density Eh depends on the shape of 
the density PDF of the I SM, which may well depend on the 
gas surface density (e.g., lKravtsovll2003h . This dependence 
may not be captured by a single clumping factor independent 
of Eh. 

The strong trend of the KS relation with metallicity has 
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Fig. 6. — Surface brightness profiles of the most massive galaxy at z = 2.5 in 
the self-consistent simulation (i.e., the simulation for which the KS relations 
shown in Fig.|3] red solid line), and in the test run in which dust abundance 
and FUV flux were set to their Milky Way values at z = 4 and the simulation 
was run with these fixed values to z = 2.5 (KS relation for this run is shown 
in Fig.|4] dashed blue line). The bottom A:-axis is in physical kpc, while the 
top axis shows the corresponding angular scale in arcseconds. The surface 
brightness, /ipeoew. w^s computed using the Flexible SPS code of Conroy & 
Gunn (2009) in AB magnitudes for the ACS F606W filter by appropriately 
redshifting the spectrum of emission to z = 2.5 (the light in this filter is thus 
dominated by the rest-frame FUV radiation) and takes into account cosmo- 
logical (1 + 2)* dimming, but does not take into account any dust absorption. 
For comparison, the dotted lines show exponential disk profiles with Re = I 
and 2 kpc. 

several interesting implications for theoretical scenarios of 
galaxy formation and interpretation of observations of high 
redshift galaxies The trend implies, for example, that star 
formation in lower-metallicity high redshift galaxies should 
be more centrally concentrated compared to their local coun- 
terparts. In particular, the regions of high redshift disks with 
Eh < 100 MqPC"^ should be considerably dimmer than ex- 
pected from the z = KS relation. 

This trend is illustrated in Figure s|5] and |6] which show dis- 
tribution of young stars and rest-frame FUV surface bright- 
ness profiles (c omputed using the Flexible SPS code of 
IConroy & Gunnl2009 ) in the most massive galaxy in our self- 
consistent simulations compared to the test run, in which dust- 
to-gas ratio and UV flux were fixed to their MW values be- 
tween z = 4 and z = 2.5. The galaxies in this test run exhibit KS 
relation consistent with the local measurements (see Fig. Hji. 
The star formation in the self-consistent run is concentrated 
towards the higher gas surface density central regions. 

Figure |6] shows that the apparent size of the simulated 
galaxy sensitively depends on the physics of molecular gas 
and the form of the KS relation. In both self-consistent and 

For example, significantly different KS relation for high redshift galax- 
ies may strongly bias gas mass estimates that use z = calibration of this re- 
lation for systems of metallicity significantly lower than solar (e.g.. lErb et all 
i2006b;.Mannucci et al...2009a) . 
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test runs, the rest-frame UV emission of the simulated galaxy 
is well described by the extended exponential disk profile 
(dotted lines). The effective half-light radii of the disks, how- 
ever, are different by a factor of two: = I kpc and 2 kpc in 
the self-consistent and test runs, respectively. This illustrates 
that to predict the observed sizes of galaxies, it is not suffi- 
cient to model dynamics of gas condensation and acquisition 
of angular momentum correctly. Correct modeling of star for- 
mation in gas of different densities is equally important. 

Incidentally, the search for extended low surface brightness 
emission in the HUDF carried out by Wolfe & Chen (2006) 
did not identify any objects that had isophotal radii /?28.4 > 
4 kpc corresponding to /iF606w = 28.4. Figure |6] shows that 
surface brightness profile of the galaxy in the self-consistent 
simulations is consistent with the observed paucity of such 
extended objects in the HUDF, while the profile of the galaxy 
in the test run is not consistent with that search. This is an- 
other illustration of the fact that the local KS relation would 
predict existence of extended low surface brightness disks at 
z ^ 3 that are not ob served in the HUDF, as emphasized by 
lWolfe&Chenl(l2006l) . 

The emerging picture based on observational constraints 
JWolfe & Chen 2006; Baker et al. 2004b; Wolfe et al. 2001 
12008) and our simulations is therefore that DLA systems cor- 
respond to the lower metallicity low-mass halos in which 
there is little or no ongoing star formation or to the outskirts 
of low metallicity disks surrounding centrally-concentrated, 
actively star forming Lyman Break Galaxies. The latter pic- 
ture is qualitatively similar to the low-metallicity z = dwarf 
galaxies, which often have a very compac t stellar system em- 
bedded in an extend ed gas-rich disk (e.g.. lMeurer et al.lll996l ; 
iBegum et al.l 1200 ^. In this picture the neutral atomic gas 
probed in absorption is largely inert to star formation, but ac- 
tively star forming regions providing metals and UV heating 
may be located in the center or in spiral arms a few kiloparsecs 
away within the same system. 

The long gas consumption time scales at Sh 100 Mq pc"^ 
would allow newly accreted gas to accumulate at these rela- 
tively low gas surface densities, instead of being consumed 
on a « 1 - 2 Gyr time scale. This can result in a more ex- 
tended star formation history if the host galaxy evolves quies- 
cently or provide gas reservoir for a starburst if the host object 
undergoes a major merger It is also possible that accumu- 
lating gas in the disk undergoes a starburst as it is enriched 
by in situ star formation and its surface density increases 
due to accretion of new gas, thereby creating conditions for 
gravitational instabilities and wides pread star formation (e.g., 
iBournaud etal]l2007l; lAgertz et al.l 2009; Dekel et al. 2009). 
The overall effect of considerably suppressed star formation 
at Eh ^100 M0 pc~^ is to preserve more gas in smaller mass 
objects for star formation in galaxies forming the bulk of their 
stars at lower (z < 2) redshifts. 

Another potentially very important implication of the 
more centrally concentrated star formation is that the stars 
formed at high-redshifts are confined to the high surface 
density regions and are therefore more resistant against 
dynamical heating in mergers. The outer gaseous disks 
are resilient to major mergers (e .g.. [Robertson et al.l |2004 
I2OO6I; ISprin gel & H ernquisll 12005]) and would help maintain 
nascent stel lar disks dynamically cold during minor mergers 
jMoster et al.i i2009ai) . thereby facilitating formation of thin 
stellar disks at later epochs and helping explain the prevalence 
of such disks in the local universe. 
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Fig. 7. — The stellar mass formed in the last 500 Myr in the simulated 
galaxies as a function of the total virial mass of their host halo for the full 
cosmological run (red circles) and the test run (blue squares). To get a rep- 
resentative sample, we show thi'ee simulation outputs together, at a = 0.24, 
0.26, and 0.28. These outputs are separated by about 250 Myr and, thus, only 
partially correlated. 

All these effects should be particularly salient in dwarf 
galaxies (Mtot ^ 10" Mq), which have the lowest metallic- 
ities and surface gas densities and where suppression of star 
formation could thus expected to be particularly strong. We 
do indeed see a progressively lower efficiency of star forma- 
tion in halos of lower mass (Fig. [U, and this decrease is not 
due to supernova feedback (as demonstrated by a comparison 
with test simulations in which the supernova feedback is arti- 
ficially disabled). The fact that stellar mass is significantly de- 
creased for a given Mtot in our star formation model compared 
to the standard star formation prescription is demonstrated 
in Figure |7] This figure shows masses of stars in simulated 
galaxies younger than 500 Myrs versus the total virial mass 
of their host halos in the full self-consistent simulation and in 
the test run in which dust-to-gas ration and FUV field were 
fixed to their Milky Way values. The latter test is equivalent 
to simply using the fixed density threshold for star forming 
gas, as is usually assumed in galaxy formation simulations. 

The figure shows that self-consistent simulation with 
our H2-based star formation model produces systematically 
smaller stellar masses in simulated galaxies compared to the 
test run. The difference is only « 20 - 30% for the largest 
galaxies, but increases to a factor of two and more for smaller 
galaxies. Remarkably, for halo masses of Mtot ^ lO'^Mo the 
star formation is completely suppressed in the full simulation, 
while in the test simulation star formation is unabated down 
to smallest resolved halo masses, well below lO'^M©. This 
shows that low dust abundance and high FUV field can play 
a major role in suppressing star formation in low metallicity, 
low mass halos - the effect that can help explaining the dif- 
ference between predicted steep small-mass slope of the halo 
mass function and the shallow observed faint-end slope of the 
luminosity function. Note that this suppression is achieved 
without efficient supernova feedback. 

The results we presented in this paper thus indicate that 
other processes in the ISM in addition (or perhaps in lieu) 
of the standard supernova feedback may be shaping the evo- 
lution of galaxy population and setting the observed disk 
sizes. Investigating these processes is an interesting avenue 
for future research. Theoretical investigations in this direction 
would be greatly aided by the observations of molecular gas in 
high-redshift star forming galaxies provided by the CARMA 
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and ALMA observatories in the near future. 
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